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Abstract 

In this paper, wc consider a contact problem with adhesion between a viscoelastic body and 
a rigid support, taking thermal effects into account. The PDE system we deal with is derived 
within the modelling approach proposed by M. Fremond and, in particular, includes the 
entropy balance equations, describing the evolution of the temperatures of the body and of 
the adhesive material. Our main result consists in showing the existence of global in time 
solutions (to a suitable variational formulation) of the related initial and boundary value 
problem. 

Key words: contact, adhesion, entropy balance, thermoviscoelasticity, global in time existence 
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1 Introduction 

This paper addresses the analysis of adhesive contact between a viscoelastic body and a rigid 
support, in the case when thermal effects are included. Contact with adhesion is described 
using the modelling approach proposed by Fremond (see [211 Chap. 14]), which was originally 
introduced for the isothermal case, combining the theory of damage (see, e.g., [22], [12], [2T] 
Chap. 12]) with the theory of unilateral contact. Indeed, although the unilateral contact theory 
(which prescribes the impenetrability condition between the bodies) does not allow for any 
resistance to tension, in the adhesion phenomenon resistance to tension is given by micro-bonds 
on the contact surface, preventing separation. Adhesion is active if these bonds (one may think 
of a "glue" on the contact surface) are not damaged. Thus, the description of this phenomenon 
has to take into account the state of the adhesive bonds (through a "damage parameter") and the 
microscopic movements breaking them, as well as macroscopic deformations and displacements. 

In the recent papers [1] and [5], we have introduced the model and derived the corresponding 
initial and boundary value problem in the isothermal case. The resulting PDE system couples 
an equation for macroscopic deformations of the body and a "boundary" equation on the contact 
surface, describing the evolution of the state of the glue by a surface damage parameter. The 
system is highly nonlinear, mainly due to the presence of nonlinear boundary conditions and 
nonsmooth constraints on the physical variables. In [3], existence of a global in time solution for 
a weak version of the corresponding PDE system was proved in the case of irreversible damage 
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dynamics for the glue. In the subsequent contribution [5], focusing on the reversible case, we 
proved well-posedness results and further investigated the long-time behaviour of the solutions. 

In this paper, we aim to generalize the model introduced in [1] and [3, including thermal 
effects both on the contact surface, and in the interior. We believe this to be interesting from 
the modeling perspective, because external thermal actions can in fact influence the state of the 
adhesive material, see |21j . 

In extending the model to the non-isothermal case, we shall adopt the following viewpoint: 
we shall assume that the body temperature and the glue temperature may be different and thus 
governed by two distinct entropy balance laws. 

1.1 The model and the PDE system 

Let us now introduce the model and derive the corresponding initial and boundary value problem. 
On a time interval (0, T), we investigate the mechanical evolution of a thermoviscoelastic body 
located in a smooth bounded domain C M^, whose boundary is 50 = Fi U r2 U Tc- Here Fj, 
i = 1, 2, c, are open subsets in the relative topology of 90, each of them with a smooth boundary 
and disjoint one from each other. In particular, Tc is the contact surface. Hereafter we shall 
suppose that both and Fi have positive measure. Without loss of generality , we shall treat 
Fc as a flat surface and identify it with a subset of . 

The thermomechanical equilibrium of the system is described by the state variables. We 
consider the absolute temperature of the body and the symmetric linearized strain tensor e(u) 
(we shall denote by u the vector of small displacements), defined in O x (0,T). Moreover, the 
variables describing the equilibrium on Fc x (0, T) are the absolute temperature of the glue 
a damage parameter X, its gradient VX, and the trace U|j, of the displacement u on the contact 
surface. The parameter X is assumed to take values in [0, 1], with X = for completely damaged 
bonds, X = 1 for undamaged bonds, and X G (0, 1) for partially damaged bonds. 

The free energy of the system is given by a volume contribution and a surface one ^Vc- It 
is known from thermodynamics that the free energy is concave with respect to the temperature. 
Thus, considering a fairly general expression for the purely thermal contribution in the free 
energy (cf. [Tl]) and normalizing some physical constants, we assume in x (0, T) 

"fn = -J W + p{^Me{u)) + ie(u)Ke(u), (1.1) 

where j is a sufficiently regular, increasing, and convex real function, the function p accounts 
for the thermal expansion energy, and K = {aij^h) denotes the elasticity tensor for a possibly 
anisotropic and inhomogeneous material. In the sequel, for the sake of simplicity, we assume 
p{'&) = "i?. A few comments on the function j are now in order. A possible choice for j, often 
used in the literature, is 

j(^) log (1.2) 

This enforces the physical constraint that t!) be strictly positive. However, in our mathematical 
analysis we are going to tackle more general situations. In particular, we shall not require any 
condition on the domain of j. Analogously, we prescribe in Fc x (0,T) 

^r. = -j(^s) + A(X)(^, - T9eg) + + a{X) + ^\VX\^ + ^X|u|^^ ^ + /_(u|^^ • n), (1.3) 

where "deq > is a critical temperature and A is a regular (quadratic) function. Once we consider 
contact with adhesion as the effect of a phase transition between the undamaged and damaged 
state of the adhesive substance on the contact surface, A' formally corresponds to the so-called 
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latent heat in phase transitions models and i^eq to the critical temperature between undamaged 
and damaged adhesion. This relates to the assumption that just by temperature devices we 
can damage the micro-bonds on the contact surface. Moreover, the indicator function I[o,i] of 
the interval [0,1] accounts for physical constraints on X, being I[o,i](^) = if X G [0,1] and 
-^[0,1] (-^) = otherwise. Analogously, denoting by /_ the indicator function of the interval 
(—00,0], the term /_(u|^ • n) renders the impenetrability condition on the contact surface, as 
it enforces that U|j, • n < (n is the outward unit normal vector to Tc). Finally, the function 
a is sufficiently smooth and possibly nonconvex, being related to nonmonotone dynamics for X 
(from a physical point of view, it corresponds to some cohesion in the material). 

The free energy describes the thermomechanical equilibrium of the system in terms of fixed 
state variables. Hence, we follow the approach proposed by J.J. MOREAU to prescribe the dissi- 
pated energy by means of a dissipation functional, the so-called pseudo-potential of dissipation, 
which is a convex, nonnegative functional, attaining its minimum when the dissipation (de- 
scribed by the dissipative variables) is zero. The dissipative variables defined in 17 x (0, T) are 
V"!? and e(ut). Thus, we define the volume part <I>q of the pseudo-potential of dissipation by 

^n = lm' + le{ut)K,e{nt), (1.4) 

where Ky = (bijkh) denotes the viscosity tensor for a possibly anisotropic and inhomogeneous 
material. The surface part of the pseudo-potential of dissipation depends on V'&s, ^t, and 
also on the difference (i9|j, — i^s) between the temperatures of the body and of the glue on the 
contact surface, namely 

<^r. = ^|V^,p + ^|Xi|2 + h{X){^\^^ - ^,)2. (1.5) 

Here, /c is a sufficiently regular function and its physical meaning is related to the heat exchange 
between the body and the adhesive material. It is of fairly natural evidence (see also [151 SOI [31] 
that the possibility (and the amount) of heat exchange between the body and the contact surface 
depends on the fact that the adhesion is more or less active. We let k to be nonnegative (in 
accordance with thermodynamical consistency ensured by the convexity of the pseudo-potential 
of dissipation), increasing, and possibly vanishing when X attains its minimum value 0. Indeed, 
we may think that if the adhesion is not active no heat exchange is allowed {k(0) = 0) or 
that a residual heat exchange is preserved even for the completely damaged adhesive substance 
(M0)>0). 

Now, let us introduce the equations in accordance with the laws of thermomechanics. We 
consider the momentum balance (in the quasi-static case) 

-divi; = f inl7x(0,r), (1.6) 

where S is the stress tensor, combined with the boundary conditions (R is the reaction on the 
contact surface) 

Sn = R inrcx(0,r), u = inrix(0,T), Sn = g inr2x(0,T), (1.7) 

f being a volume force and g a traction. The thermal balance is given by the following entropy 
equation 

st + dwQ = h inOx(0,r), (1.8) 

s denoting the internal entropy, Q the entropy flux, and h and external entropy source. Indeed, 
equation (jl.Sp can be obtained rescaling the first law of thermodynamics (dividing the internal 
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energy balance by the absolute temperature), under the small perturbation assumption. We 
refer to [TU], [TT] and 0, [H]) [S] for details on this modelling approach and related analytical 
results. Moreover, we supplement (jl.Sp with the following boundary conditions 

Q -n = F on Tc X (0,r), Q • n = on \ x (0, T), (1.9) 

where F denotes the entropy flux through Tc- Hence, on the contact surface, we introduce a 
balance equation for the microscopic movements, that is 

S-divH = inTc X (0,r), H • n, = on x (0, T), (1.10) 

B, H representing interior forces, responsible for the damage of adhesive bonds between the body 
and the support, and the outward unit normal vector to dTc- Then, the entropy equation on 
the contact surface is given by 

(9tSs + divQ^ = F in Tc x (0,r), • = on x (0, T). (1.11) 

Here, Sg is the entropy on the contact surface, Qs the surface entropy flux, and the term F, 
given by the flux through the boundary Tc (cf. p.9|) ). represents a surface entropy source. 

Constitutive relations for S, R, s, Q, F, H, s^, Qs are given in terms of the free energies 
and the pseudo-potentials of dissipation. More precisely, we have 

^ ^(??)-divu, (1.12) 



^=1{{^,)-X{X), (1.13) 



where ^ is the derivative of the convex function j. In the physical case j{x) = xlogx — x (cf. 
()1.2p ). we have 

i{x)=\ogx. (1.14) 

In particular, I in (11.140 in fact yields an internal positivity constraint on the system tempera- 
tures "d and "dg- Furthermore, 

«'=-|ls: = -^''" 

The constitutive relation for the stress tensor S accounts for dissipative (viscous) dynamics for 
deformations, in that we have 

S = 777^ + Try^ = Ke{n) + K„e(ui) + ^1, (1.18) 
(9e(u) de{\it) 



(1 denotes the identity matrix), while the reaction R is given by 

dvL\ 



R = = -Xu|^^ - a/_(u|^^ . n)n. (1.19) 
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We further prescribe B 



d^fr d<^r 



2' 



{dl- and i9/[o,i] standing for the sub differentials of I_ and -/^[o,i]) respectively), and let H be 

H = — -^ = VX. (1.21) 

Remark 1.1. Let us point out that the evolution of the system is characterized by dissipation 
due to choice of dissipative potentials <I>q and <I>rc (cf- (|1-4|) and (|1.5|) ) and the balance laws of 
thermodynamics (see in particular (jl.Sp and (jl.lip . ()1.6p . and (ll.lOh ). 

Combining the previous constitutive relations with the balance laws, we obtain the following 
boundary value problem 



dtiii'd)) - div(ut) -A^ = h in U x (0, T), 

'o in (ao\rc) X (o,r), 

^-k{X){^-^s) inr, x(o,r), 

dMi^s)) - dt{X{x)) - Ai)s = Kx){^ - inr^x (o,r), 

dn^s = inar, x(0,T), 

- div {Ke{u) + K^e{ut) + ^1) = i in x (0, T) , 

u = inrix(0,T), {Ke{u) + K^e{ut) + 'dl)n = g inr2x(0,T), 
{Ke{u) + K^e{ut) + ^l)n + Xu + dI-{u -11)113 in Tc x (0,r), 

Xt - AX + dI[o,i]{X) + a'{X) - A'(X)i?eq ^ -A'(X)i?, - ^jup in T, x (0,r), 
dnX = in dTc x (0, T) 



1.22) 

1.23) 

1.24) 
1.25) 
1.26) 
1.27) 
1.28) 

1.29) 

1.30) 



(here and in what follows, we shall omit for simplicity the index v^^ to denote the trace on Tc 
of a function v, defined in fi). 

Remark 1.2. As we have already pointed out, our analysis actually accounts for a form of the 
thermal contribution in the free energies more general than (ll.2h (cf. (jl.lh and (11.30 ). Thus, 
our results may apply to several physical situations with different thermal behaviour. We note 
that we can also handle the case when the specific heat cy, given by the thermodynamic relation 

cv = (cf. Illl Rem. 2.11), is not constant. Indeed, it is known from physics that the 

specific heat may be depending on the temperature, e.g., cv(i?) = i?'^, 7 > 0. A choice of j 
corresponding to the latter expression of cy is 

^7+1 

m 



7(7 + 1) ' 



which is covered by our analysis. In particular, letting 7 = 1 (and hence j{'d) = 2")-, we get 
£(??) = 19, so that p.22p and (|1.24p reduce to Caginalp-type heat equations (cf. with |16j). 
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1.2 Related literature and our own results 

For a review of the theory of contact problems see, e.g., the monographs [271 [El EH], and the 
references therein. We also refer to [H [5] for some partial survey of the literature on (isothermal) 
models of adhesive contact. In this connection, we mention [21], where a unilateral contact model 
(also derived within Premond's approach) is considered. Therein, the adhesive properties are 
described by the condition Xu = on the contact surface. The author proves an existence 
theorem for the related PDE system and also develops some numerical investigations. The 
paper |25j focuses on a model combining unilateral contact with adhesion and friction as well: 
under a smallness condition on the friction coefficient, an existence result is proved and various 
numerical schemes are proposed. 

As for the literature on contact models including thermal effects in the three-dimensional 
framework, besides the contributions mentioned in [TH] we may recall where a frictionless 
contact problem between a thermoelastic body and a rigid foundation is modelled by a parabolic 
equation for t?, coupled with an elliptic equation for u, with mixed boundary conditions. An exis- 
tence result is proved, provided the coefficient of thermal expansion is sufficiently small. Among 
dynamic models, in which contact is rendered by means of a normal compliance condition, we 
quote [20] , which deals with a wide class of frictional contact problems in thermoelasticity and 
thermoviscoelasticity. Moreover, we mention [1], where a frictional contact problem involving 
a thermoelastic body undergoing wear on the contact surface is investigated. A well-posedness 
result is proved for a system coupling a parabolic equation for the temperature, a variational 
inequality for the displacement, and a first order equation for the wear function, supplemented 
with nonlinear boundary conditions. 

Nevertheless, as far as we know, no results are available, in the literature, on unilateral 
contact models which take into account both adhesive properties and thermal effects. Indeed, 
one of the main novelties of the present contribution is that we consider heat generation effects 
in the adhesive contact phenomenon, too. That is to say, we allow for the body and the adhesive 
material to have different temperatures, whose evolution is mainly ruled by the heat exchange 
throughout the contact surface. More precisely, the entropy flux F through Tc (occurring in 
()1.9p ) plays the role of a source of entropy in (II. lip . From an analytical point of view, this 
results in a nonlinear coupling between (jl.22p - ()1.23p and (jl.24p and gives raise to some technical 
difficulties. A further peculiarity of our work consists in assuming entropy balance laws (in place 
of the more usual internal energy balance), for describing the evolution of the body and of the 
glue temperatures. This brings to strong and possibly singular nonlinearities in ()1.22p and (ll.24p 
(see (I1.14P ). An advantage of this choice is that, assuming that the domain D{j) C (0, +oo) 
(as in the case of the classical choice (|1.2p ) once the problem is solved in a suitable sense, the 
positivity of the temperatures is deduced. On this fact the thermodynamical consistency of the 
model relies, see also [1^ [71 [SI [TT]. This is of particular interest in the present case, since the 
low spatial regularity of the solution components i? and i?^, along with the nonlinear boundary 
condition (jl.23p . prevents us from using any maximum principle. 

In fact, we shall study the Cauchy problem for a generalized version of system (ll.22p - ()1.30p . 
see Problem (P) in Sec. 12.31 Namely, we replace the subdifferential operators in (|1.28p - (|1.29p . by 
general maximal monotone operators (possibly rendering physical constraints on the variables X 
and u). Further, we generalize the choice of the nonlinearity i, allowing for a maximal monotone 
operator. In fact, the only restriction we impose on £ is that the resulting internal energy of 
the system be coercive, cf. with ()1.3ip below. This is reasonable from a physical point of view 
and still enables us to include several choices of i in our analysis (in particular, (|1.2p . as well 
as the examples of Remark ll.2p . The idea is that, once the internal energy of the system is 
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bounded, the absolute temperature is bounded, too. Let us focus on the volume temperature 
i9. As known from thermodynamics, the internal energy (depending on the entropy s) can be 
introduced as the convex conjugate function (with respect to the variable •&), of the negative of 
the free energy (which is convex w.r.t. the temperature as the function j is convex, see 
and ()1.3p ). Namely, in the case of the volume free energy the related internal energy is 

e{s, ■) = i-^nYis, •) = snp{s^ + ^n{^, ■)). 

Thus, our coercivity condition may be expressed in terms of the conjugate of j by 

3Ci,C2>0: f{y)>Ci\x\-C2 if y = £{x) (1.31) 

(see subsequent ()2.H2p and Remark 12. 2p . 

The main difficulties attached to the analysis of the PDE system (jl.22p - ()1.30p are related 
to the singular character of the entropy equations (|1.22p and p.24p . to the nonlinear coupling 
between the latter equations, as well as between (I1.28P and (ll.29p . and, last but not least, 
to the presence of general multivalued operators in all of the latter equations. In particular, 
it seems to us that dealing with a general maximal monotone operator in (ll.22p and (ll.24p 
brings about some technical difficulties, particularly in connection with the third type boundary 
condition (I1.23P for -d on Tc- 

All of these peculiarities will be carefully handled in the proof of our main result. Theorem 
1 (see Sec. 12. 3p . stating the existence of global in time solutions to the Cauchy problem for (the 
generalized version of) system (jl.22p - (jl.30p . We sketch below the main steps of our procedure, 
based on a suitable approximation of Problem (P), and on the derivation of suitable a priori 
estimates, which enable us to pass to the limit in the approximation. Such estimates are intrin- 
sically related to the dissipative character of the system, highlighted in Remark 1 1.1[ In fact, in 
the paper [6] we take advantage of the dissipative character of the system to perform its long- 
time analysis, showing that in the limit we reach a stationary equilibrium in which dissipation 
vanishes. 

Note, however, that uniqueness is still an open issue, at least in the functional framework of 
our existence theorem. Without going into details, we may point out that the major obstacle 
is due to the singular character of equations (jl.22p and (|1.24p . In particular, the boundary 
condition (ll.23|) makes it harder to prove contraction estimates leading to uniqueness. We refer 
to Remark l3.2l for additional observations on this point. Actually, uniqueness holds in the (more 
regular) framework of the approximate problem, see Section [3.51 

Plan of the paper. In Section 2, we enlist our assumptions on the problem data, present the 
variational formulation of the Cauchy problem for (a generalized version of) system (jl.22p - (jl.30p . 
and state of our main result. In Section 3, we set up the approximation of Problem (P), suitably 
regularizing the maximal monotone operators in equations (jl.22p and (|1.24p and therein insert- 
ing (vanishing) viscosity terms. Hence, we prove a well-posedness result for the approximate 
problem. We combine a Schauder fixed point technique for local existence with a prolongation 
argument, based on global in time a priori estimates, while uniqueness follows from contraction 
estimates. Next, in Section 4 we pass to the limit in the approximate problem by compact- 
ness and monotonicity tools, and show that the approximate solutions converge to a solution of 
Problem (P). Finally, in the Appendix we prove some auxiliary technical results. 
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2 Main result 

2.1 Setup and preliminary results 

Notation 2.1. Throughout the paper, given a Banach space X, we denote by x'{'i')x the 
duality pairing between X' and X itself, and by || • ||x both the norm in X and in any power 
of X; C^{[0,T]; X) is the space of weakly continuous X-valued functions on [0, T]. Whenever 
X = Yi X . . . X Yat, we denote by vTj, i = 1, . . . , the projection on the i-th component. 

Young inequalities. We recall the Young inequality for convolutions, namely 

Vp, g, r € [l,ool s.t - = - + --1 y a e LP(0,T) b e L^O,T; X) we have 

r p q (2.1) 

a * 6 G L''(0, T; X) and \\a * 6||l'-(o,T;X) < \\a\\LP(o,T) \\b\\Liio,T;X) , 
and the Young inequality 

V(5>0 3C5>0: Vp, g G (l,cx)) with - + - = 1 ab<6aP+Csb'' for ah a, 6 G M . (2.2) 

P Q 

Functional setup. Henceforth, we shall suppose that is a bounded smooth set of ffi'^, such 
that Tc is a smooth bounded domain of M?, and use the notation 

H := L'^in), V:=H\n), and 
W := {v G : v = a.e. on Ti} , 

the latter space endowed with the norm induced by V. We shall work with the standard Riesz 
operator 

^:V^V' given by y,{Jl{u),v)y := uv + VuVv for all n, z; G (2.3) 

Jn Jn 

and denote by IRfc the analogously defined Riesz operator mapping H^^Tc) into {H^^Tc))'- 
Further, we shall extensively use that 

V C L^{Tc) with a continuous (compact) embedding for l<p<4 (1 < p < 4, resp.), (2.4) 
H^{Tc) C L^{Tc) with a compact embedding for 1 < p < oo. (2.5) 

For simplicity, we denote by /p uv (Jp_^ uv, resp.) the duality pairing (•j:^-i/2(r^))3 (ti, t')(j:/i/2(rc))3 
between (iJ-V2(r^))3 and {H^/^{Tc)f (between {H-^/'^{T2)f and {H^/'^{T2)f, resp.). Finally, 
given a subset C M^, iV = 1, 2, 3, we shall denote by |0| its Lebes gue measure. 

Preliminaries of viscoelasticity theory. We now introduce the standard bilinear forms of 
linear viscoelasticity which allow us to give a variational formulation of equation p.26p . Dealing 
with an anisotropic and inhomogeneous material, we assume that the fourth-order tensors K = 
(o-ijkh) and Ky = {bijkh}-, denoting the elasticity and the viscosity tensor, respectively, satisfy 
the classical symmetry and ellipticity conditions 

O-ijkh — (^jikh — Oikhij i bijkh — bjn^ji — b^fiij , i,jj k,h — 1, 2, 3 

3 ao > : aijkhiijikh > ao(.ij(.ij ■ = ^ji , ^, i = 1, 2, 3 , 

3(3o > : bijkhiijikh > f^oiijiij '^s.ij ■ = iji , j = i, 2, 3 , 
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where the usual summation convention is used. Moreover, we require 

aijkh,bijkh G L°°{i^) , i,j,k,h = 1,2,3. 

By the previous assumptions on the elasticity and viscosity coefficients, the following bilinear 
forms a, b : W x W — > M, defined by 

a(u,v):= / aijfc/iefc/j(u)ejj(v) Vu,veW, 
Jn 

6(u,v):= / bykh^khM^iji^) Vu,vGW 
Jn 

turn out to be continuous and symmetric. In particular, we have 

3M > : |a(u,v)| + |6(u,v)| < M||u||wi|v||w Vu,vGW. (2.6) 

Moreover, since Fi has positive measure, by Korn's inequality we deduce that a(-, •) and •) 
are W-elliptic, i.e., there exist Ca,Ch > such that 

a(u,u) > a||u||^ VuGW, (2.7) 
6(u,u) > a||u||^ VuGW. (2.8) 

Relying on Green's formula (see, e.g., [18]), the variational formulation of (11.260 (cf. (12.360 
below) can be derived by a standard argument. 

2.2 Statement of the assumptions 

As we mentioned in the introduction, we shall address the analysis of a generalized version of 
system (ll.22p - (ll.30p . in which the operators occurring in (11.270 and (11.290 are replaced by gene- 
ral maximal monotone operators. We now enlist our assumptions on the involved nonlinearities 
and on the problem data. 
We consider 

a proper, convex, and l.s.c. function j -.M ^ {—oo, +oo] , (2-Hl) 

and 

its sub differential in the sense of convex analysis i = dj : M ^ 2* . (2.9) 
The crucial assumption on j is that the following coercivity condition holds: 

3Ci, C72>0 yx eD{i), y £e{x) : yx - j{x) > Ci\x\ - C2 . (2.H2) 

Remark 2.2. Note that (j2Jl2l) can be rephrased as 

3Ci, C2>0 Vx G D(£), y G ^(x) : j* (y) > Ci\x\ - C2 . (2.10) 

In particular, (|2.10p is fulfilled when j{x) = x{log{x) — 1) for all x G (0, +00) and hence i is the 
logarithmic nonlinearity, i.e. £{x) = log(x) for all x G (0, +00). In this case, simple computations 
show that j*{y) = e?^ = i~^iy) for all y G M, whence (j2J0|) . 

Henceforth, we shall denote by 

7 the inverse of the operator i (recall that 7 = dj*) 

and, with a slight abuse of notation, we shall call £ as well the realization of (12.90 as a 
maximal monotone operator i : L^{0,T; H) — > 2^ io,T;H) ^ maximal monotone operator 
i : L2(0,T;L2(rc)) ^ 2^'(0'^'-^'(r,))^ respectively). 
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Assumptions on the other problem nonhnearities. Further, we let 

a : {H^^'^{rc))'^ — > [0, +00] be a proper, convex and l.s.c. functional, 

with 3(0) = = minS, and (2.H3) 
we set a:=da: {H'/\r,)f ^ 2(^"''(r.))3 ^ 
Indeed, a shall generalize the subdifferential operator appearing in (11.280 . 

Remark 2.3. In fact, condition (jl.28p may be rendered rigorously (see [H Sec. II] for details) 
by introducing the set 

X„ := {v G (i/i/2(rc))3 : V n < a.e. in TJ, 

with indicator function /x_ . Then we consider the maximal monotone operator 9/x_ : {H^^^{Tc))^ 
2(^-'''(^^^^y\ given hy 

T] G {H-^/^{rc)f belongs to dlx_{y) if and only if 
yGX_, / r7-(v-y)<0 Vv G X_. 



In the same way, in equation (jl.29p we shall consider 

a maximal monotone operator /3 : M ^ 2^*, with domain D(/3) C [0,+cx3). (2. 114) 



It is well known that there exists a proper, l.s.c. and convex function fd : D(/3) (—00, +00] 
such that P = dp. Concerning the functions A, k, and a', we assume that 

o"' : M — > M is Lipschitz continuous, with Lipschitz constant L^^, (2. 115) 

k : M — > [0, +00) is Lipschitz continuous , with Lipschitz constant L^, (2. KG) 

AgC^'^(M), (2.H7) 

i.e., A has a Lipschitz continuous derivative. As a consequence, 

3C3 >0 Vx,2/GM : |A(x)-A(y)| <C3(|x| + |y| + l)|x-y|, (2.11) 

3C73>0VxGM: |A'(x)| < C^dxj + 1) . (2.12) 

Assumptions on the problem data. We assume that 

heL^ (0, T;V')n (0, T; H) , (2.H8) 

i e L^{0,T;H^), (2.H9) 

g£L\0,T;{H~^/\r2)f). (2.H10) 
Then, we remark that the function F : (0, T) — > W' defined by 

w/(F(t),v)w : / f(t) •v+ / g{t)-v VvGW^ for a.e. t G (0, T), 

fulfils 

F G L2(0,T;W'). (2.13) 
Finally, we require that the initial data fulfil 

woeH with j*{wo) G L^{n) , (2.14) 

zo G L\T,) and f{zo) G ^^(r,) , (2.15) 

uo G W and uq G D(S) , (2.16) 

XoeHHr,), p{Xo)£L\r,). (2.17) 
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2.3 Variational formulation and statement of the main result 

We now state the variational formulation of the initial-boundary value problem for a generalized 
version of system (jl.22p - (|1.30p . featuring the nonlinearities introduced above. 

Problem (P). Given a quadruple of initial data {wq, zqi uq, Xq) complying with (|2.14|) - (j2.17|) . 



find w, "ds, -2) u, X, r], ^) such that 

1? G L2(0,r;y) nL°°(0,T;L^(17)), (2.18) 

w e L^{0,T;H)r\H^{0,T-V'), (2.19) 

f{w)eL°^iO,T;L\n)), (2.20) 

G L\0, T- H\T,)) n L°°(0, T; L\T,)) , (2.21) 

z e L°°{0,T;L\r,))nH\0,T;H\T,y), (2.22) 

j*iz)eL'^iO,T;L\T,)), (2.23) 

ue H\0,T;W), (2.24) 

r/eL2(0,r;(i/"V2(p^))3)^ ^2.25) 

X e L2(o,r;/?2(r^))nL~(o,r;//i(re))ni7i(o,T;L2(re)), (2.26) 

^eL\0,T;L\T,)), (2.27) 
fulfilling the initial conditions 

w{0)=wo a.e. inn, (2.28) 

z{0) = zq a.e. in , (2.29) 

X(0) = Xo a.e. in , (2.30) 

u(0) = uo a.e. in (2.31) 

and 

y,{Wt,v)y- I div(Ui)w+ / \/^Vv+ I k{X){'d-^s)v 

Jn JU JTc (2.32) 

= Yi{h,v)y yv£V a.e. in (0, T) , 

w{x,t) e £{'d{x,t)) for a.e. (x, t) e X (0, T) , (2.33) 

Jr, JTc ^2.34) 
= / k{X){i} - i}s)v yveH^{Tc) a.e. in(0,r), 

z{x, t) G e{^s{x, t)) for a.e. {x, t) eV^x (0, T) , (2.35) 

6(ut, v) + a(u, v) + / i?div(v)+ / (Xu + 77) • v 

Jn Jfc (2.36) 

= w'(F,v)w VveW a.e. in (0,r), 

a.e. in (0, T), (2.37) 

Xt- AX + ^ + a'{X) = -X'{X)'3s-^\u\'^ a.e. in Tc x (0, T), (2.38) 

^ e /3(X) a.e. in x (0,r), (2.39) 

dn,X = a.e. in x (0, T) . (2.40) 
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Note that, to simplify notation, we have incorporated the contribution — A'(X)i?eq occurring 
in (fT^ into the term cr'(X) in (OSI) . 

Theorem 1 (Existence of a global solution). Assume (l2.Hl|) -( )2.H10p . Then, Problem (P) 
admits at least a global solution ("i?, tt), i?^, z, u, X, 77, ^) with the regularity (|2.18p - (|2.27p . 

Remark 2.4 (Positivity of the temperature). Clearly, in the case 

D(j) C (0,+oo), 

(such as for (jl.2p ). from relations (j2.33p and (j2.35p we infer that both system temperatures ?9 
and "ds are strictly positive almost everywhere in J7 x (0,T) and in Tc x (0,T), respectively. 

Strategy of the proof of Theorem [1], We shall approximate Problem (P) by suitably 
regularizing equations p.32p and p.34p . More precisely, we shall add some viscosity terms to 
both equations, and replace the operator £ therein by a Yosida-type regularization (see ()3.ip 
below). For technical reasons (cf. with Remark 13. 2p . we shall keep the viscosity parameter 
distinct from the Yosida regularization parameter, denoting by e > the former and by > 
the latter. Hence, we shall call (Pe) the initial and boundary value problem for the resulting 
approximate system and prove that it is well-posed following this outline: first in Sections 13.21 - 
13.31 we are going to prove the existence of a local solution by a fixed point argument. Next, 
in Section 13.41 we are going to extend such a solution to the whole interval (0,T), while in 
Section [3.51 we shall obtain contraction estimates leading to uniqueness for Problem (P^). 

Finally, in order to prove Theorem [T] we shall pass to the limit in Problem (Pe ) in two 
steps. First, we shall keep n> fixed and let e \ 0: in Section ITT] we are going to show that 
the approximate solutions converge, as e \ 0, to a solution of the initial and boundary value 
problem obtained by setting e = in the approximate equations (|3.13p - (|3.14p below. Secondly, 
we shall also let ^ \ and obtain in the limit a solution to Problem (P). 

Notation 2.5. Henceforth, for the sake of notational simplicity, we shall use the same symbol 
(•, •) for the duality pairings ■)w' V" ■)v ^'^'^ H'^{Tc)'^'^ ■)//i{rc)' ^^'^^ further, denote by 

the symbols C, C' most of the (positive) constants (2-41) 

occurring in calculations and estimates. 

3 Approximation 

3.1 The approximate problem 

We approximate Problem (P) by modifying equations ()2.32p and p.34p in the following way: 

• first, we shall add to ()2.32p the regularizing viscosity term e'Jl{'dt) and to (I2.34p the viscosity 
term e^rd'^t'&s), with e > 0: this shall enable us to perform enhanced regularity estimates 
on the (approximate) equations for the temperatures "& and i?^ and ultimately to prove the 
global well-posedness of the approximate system, see also Remark \'6.2\ 

• second, both in p.32p and in ()2.34p we shall replace the operator i with its Yosida-type 
regularization 

^^■.= {fild + j^)-\ (3.1) 

where > 0, Id : M ^ M is the identity function, and 7^ : R ^ M is the (/i-)Yosida 
regularization of the inverse 7 of i, see (|3.2p below. The choice of is motivated by 
technical reasons, cf. with Remarks 13.21 and 14.21 later on. 
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For the purposes of the above approximation, we recall that the Yosida regularization 7^ : M ^ M 
is a Lipschitz continuous function, given by 

J (w) ■=-(w- pJw)) for all w €R, (3.2) 

where : M ^ M is the (/i-) resolvent operator associated with 7, defined for every S M as 
the unique solution p^{w) of the inclusion 

p^{w) -w + ii-iip^iw)) 3Q. (3.3) 

We shall also deal with the Yosida approximation j* of j*, defined for every /i > by 

i*(w) :=min| *^~^* +i*{w)\ for all u- € M. (3.4) 



j/eM 1^ 2ijl 

We recall that j* S C"'^(IR), with derivative j*' = 7^, and that it fulfils 

= f I7mHI' +i* (PmH) Vu; G M. (3.5) 
It was proved in |13l Sec. 3] that the operator is well-defined on M, monotone, and that 

M — > M is Lipschitz continuous with Lipschitz constant l/^u. (3.6) 



Approximate initial data. In order to properly state our approximate problem, depending 
on the parameters e > and /x > 0, we shall need to prescribe some initial conditions for 'd 
and 'ds ■ To this aim, in Lemma 14.11 we shall construct sequences of initial data and '&^s,£^l 
satisfying 

€,^y^ j;{^M,))^L\n), (3.7) 

€ H\T,) , j; {^^{^%,)) G L\T,) , (3.8) 
and such that there exists a constant Mq > independent of e > and > with 

{^^.{^%))\\^,^^^ + (if^«,^))||^,(r.) ^ e,p>0, (3.9) 

and a constant Mq > independent of e > but possibly depending on n > with 

^'/'lldlv^ + ^'/'lKe/.llHi(r.) < Mq'' for ah e > 0. (3.10) 

For simplicity, throughout this section we shall omit the indexes e and p, in the notation for the 
initial data. The variational formulation of the initial-boundary value problem approximating 
Problem (P) then reads: 

Problem (P^). Given a quadruple of initial data (??°, uq, Xq) fulfilling (IT71) - (ITTU]) and 
(f2J6]) - (f2T71) . find functions (??, u, X, ^, r?) satisfying ([221]) -(1122]), with i} and 'ds such that 

i^eH^{0,T;V), (3.11) 
Gi7H0,T;i7i(r,)), (3.12) 
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e [ ^tv+ [ dt^f,{^)v- [ div{ut)v + e [ V^t^v + [ Vi9Vu 
Jn Jn Jn Jn Jn 

+ / k{X){'d - ^,)v = {h,v) yv€V a.e. in (0,r), 

- / dt\{X)v+ I V^s^v= I k{X){'d-^s)v yveH^iTc) a.e. in (0,r), 

and equations (|2.36p - (|2.40p . along with the initial conditions (|2.3ip and (|2.30p for u and X, 
respectively, and, for ■& and i9s, 

i?(0)=i?° mV, (3.15) 

MO) = ^^s inH^r,). (3.16) 

Remark 3.1. Combining (I3.1ip - ()3.12p with the fact that is Lipschitz continuous (cf. p.6p ). 
we conclude that for any solution (t?, i?^, u, X, ^, 77) there holds 

e L^{0,T;V)nH\0,T;H), e L^{0,T; H\r,))nH\0,T; L\T,)) . (3.17) 

Remark 3.2. Without going into details, a few comments on the approximate Problem (Pe) 
are in order. Besides regularizing the maximal monotone operators in (j2.32p and in (j2.34p . we 
have inserted in both equations a viscosity term in order to make each of (the Cauchy problems 
for) the approximate equations well-posed. To understand why, let us focus on equation (I3.13P 
(analogous considerations apply to equation (j3.14p l. 

Indeed, because of the nonlinear term dt{J^^{'&)) in (I3.13p . in order to prove that the (Cauchy 
problem for the) latter equation has a unique solution, one has to integrate it in time. Hence, as 
it will be clear from the calculations in Lemma [3. 5 1 the additional viscosity term e'3l{'dt) enables 
us to deal, in the integrated version of (j3.13p . with the third type boundary condition on On 
the other hand, the choice of the operator in (I3.13P (of in ()3.14p ). in place of the usual 
Yosida regularization of £, is due to technical reasons connected to the construction of sequences 
of approximate initial data fulfilling ()3.7p - ()3.10p . see Lemma l4.ll and Remark 14.21 later on. 

The ultimate reason why we keep the viscosity parameter e distinct from the Yosida pa- 
rameter in both approximate equations (|2.32p and (|2.34p is due to the fact that, in order to 
recover the L°°(0, T; f/')-regularity ()2.19p for the solution component w (the L°°(0, T; L^(rc))- 
regularity (j2.22p for the solution component z, respectively), one has to test some approximation 
of ()2.32p ( p.34p . respectively) by a function approximating w (z, resp.), and obtain some bound 
obviously independent of the approximation parameter. In the present framework, performing 
such an estimate on equation (I3.13P (on ()3.14p . resp.) with e > would not lead to estimates 
on ^fj,{'d) independent of the parameters e and /x, essentially because the term (e3?(T9j), ^^(t?)) 
( {e'JiY c{dt'& s) , s)) , resp.) cannot be dealt with by monotonicity arguments. That is why, 
in Section 14.21 we shall perform the crucial estimate leading to regularity (j2.19p and (j2.22p only 
after taking the limit in Problem (Pe ) as e \ with fi > fixed. 

The following theorem holds. 

Theorem 3.1 (Global well-posedness for Problem (Pg)). Under assumptions ()2.IIip -( l2.H10p . 

for any set of initial data ('i?'', Uq, Xq) complying with conditions (|3.7p - (|3.1Up and (|2.16p - 
(I2.17P and for alle, fi > there exists a unique global solution (t?, "i^g, u, X, ^, 77) to Problem (Pe). 



(3.13) 
(3.14) 
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As already mentioned, we shall first of all prove the existence of a local solution to Problem (Pe ) 
by means of a Schauder fixed point argument, which relies on auxiliary intermediate results on 
the existence and uniqueness of solutions for the single equations of (Pe). We shall prove such 
results in the following Section \'6.2\ and conclude the proof of local existence in Section 13.31 
Finally, in Section EH] we shall show that the local solution extends to a (unique, by Section [33|l 
global one. 

3.2 Fixed point setup 

For a fixed t S (0, T] (which shall be specified later on) and a fixed constant ii > 0, we consider 
the set 

St := \{^,^s,X) eL^{0,t;H) x L^/^{0,t;L^{T,)) x L^0(0, t; L6(r,)) : 

^ (3.18) 

||(t?, l?s,X)||^2(o,t;//)xL'V2(o,t;L3(rc))xLiO{0,t;L6(rc)) - ^\ ■ 

We shall construct an operator T (its definition is split in several steps), mapping 8^ into itself 
for a suitable time < T < T, in such a way that any fixed point of T yields a solution to 
Problem (P^) on the interval (0,r). Then, in Proposition 13.71 we shall prove that T : ^ 
admits a fixed point by the Schauder theorem. 

Notation 3.3. We shall denote by 

Mi, i = 1,2,3 a positive constant depending on R, on the problem data, 

on Mo (cf. ^Ml) and (cf. (|330]) ). 

but independent of e > 0, and of the fixed t G (0, T] 

Mf , i = 1, 2, 3 a positive constant depending on R, on the problem data, 

on Mo (cf. ([33])) and (cf. (IXTnD ). and possibly on e > 0, 

but in any case independent of the fixed t G (0, T] . 

Further, we shall keep to the notation (|2.41|) for all the constants which do not depend on the 
approximating parameter e. 

Step 1. We take (i?, t?s,X) G St and consider (the Cauchy problem) for system ()2.36p - ()2.40p . 
in which ■!? in (j2.36p is replaced by ■& and X and 'dg on the right-hand side of (j2.38p are replaced 
by X and 'dg-, respectively. 

Lemma 3.4. Assume (l2ll3]l . (I2ll4l) . (f2JT5l) . ^(Tm) . (l2ll9]K (l2.Hin[) . (f2l6l) . and (fZTTp . 

Then, there exists a constant Mi > such that for all {'&,i)s,X) G St there exists a unique 
quadruple (u, r/,X,^), with the regularity 

l|u||Hi(0,t;W) + ll^llL2(0,t;H-i/2{rc)) + ll?llL2(0,t;L2(r^)) (3.19) 

+ ll^llL2(o,t;//2(r^))nj;^cx)(o,t;//i{rc))n/i'i(o,t;L2(rc)) ^ ^-^i > 
complying with the initial conditions ()2.30p - ()2.3ip . solving the PDE system 

6(ut,v) + a(u,v)+ / (Xu + T7)-v (3.20) 
= (F,v)- / ?div(v) Vv G W a.e. in (0,t), 
Xt-AX + ^ + (T'(X) = -A'(X)t9, -^|up a.e. inTcX{Q,t), 
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and such that u and ij fulfil (|2.37p . X and ^ comply with (|2.39p - (j2.40p on (0,t). 

Let us point out that, in view of (12.H7p . ()3.18p . and the Holder inequahty, there holds 

\\^'W'^s\\L^O,t;L^r^)) < C'll^llLiO(0,t;Le(r,)) ll^^s |lL5/2(o,t.^3(r^)) < CR^ . (3.21) 

Then, also taking into account the fact that "d G L^(0, t; L^(rc)), Lemma 13.41 follows from [5l 
Thm. 1], to which we refer the reader. 
Now, we let 

Vt:={(u,X) GFi(0,t;W) x {L\0,t;H''{T,)) D L'^{0,t;H\T,)) D H\0,t;L\T,))) : 
l|u|ki(o,t;W) + ll^llL2(o,t;H2(rc))nL°°(o,t;_f/i(rc))n/fi(o,t;L2(rc)) ^ • 

(3.22) 

Thanks to Lemma [331 we may define an operator 



Ti : St ^ V. 



mapping every triple {'d,i!)s,X) € St into the pair (u, X) solving the Cauchy problem for sys- 
tem ()3.20p (we imply that with (u, X) the solution components r] and ^ satisfying p.37p and ()2.39p 
are uniquely associated). 

Step 2. We fix ('!?s,u, X) G 7r2(St) x Vt and consider (the Cauchy problem) for ()3.13p with 
data (i?s, u, X). 

Lemma 3.5. Assume (l2lIT]l . (l2Jl2l) . (l2Jl6]) . (l2Jl8| . and (l3Tp - (l3lil]l . 

Then, there exist M2,M2 > suc/i that for all (t?s,u, X) G 7r2(St) x Vt i/iere exists a unique 

ll^llL2(o,t;y)nL°=(0,t;Ll(Q)) < ^2, ll^?t II L2 (0,t;y) < ^2 ) (3.23) 

complying with the initial condition (I3.15|) and with 



e f dtv+ f dt^^{'d)v- [ dW{ut)v + e [ V^t^v + [ V^Wv 
Jn Jn Jn Jn Jn 

+ / k{X){i}-ds)v = {h,v) \fveV a.e. in (0,t). 



(3.24) 



Proof. For simplicity, throughout the proof we shall use the notation 

X{x, t) := k{X{x, t)) for a.e. (x, t) G L^ x (0, t) . (3.25) 

It follows from ([2Jl6]) . the regularity of X (cf. (IX^ ). and ([23]) that 

X G L°°(0,t;i?^(rc)), hence 3C G L°°(0,t;LP(rc)) Vl<p<oo. (3.26) 

In view of [171 Thm. 1], the Cauchy problem for (j3.24p has at least a solution i} G H^{0, t; V). 
In order to prove uniqueness, we let 'i?2 £ -^^(0, t; V) be two solutions of the Cauchy problem 
(13.151 [3^2^ . and set 'd : 'di — '!?2- We subtract the equation for 1)2 from the equation for 'di and 
integrate on (0,t), with <t <t. Thus, we get 

e [ dit)v+ [ {^^iMt))-^t.iMt)))v + e [ vd{t)Vv 
Jn Jn Jn 

+ [ {l*vd){t)Vv+ [ {l*xd){t)v = 
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for all V ^ V. Hence, we take f = ■!? and integrate in time: also using that the operator is 
monotone, with straightforward computations we obtain 

\n% +ej'^ \\vni + ^11(1 * ^m)\\l < wnLHv.) m * rnhnr.)- (3.27) 

In view of the Young inequality for convolutions (12.11) , we have 

\\il*m{s)\\LHr.)<s'/^\\XnLHo,s-,L^{r.))<s'/^\\nL^ VsG [0,t]. 

Hence, ([QT]) yields 



v<^J^ Mv + Cet\\X\\i^^o,t;HHr.)) I mhio,s;V) > (3-28) 

where depends on e and also on the embedding constants in ()2.4p - ()2.5p . By applying the 
Gronwall Lemma (see, e.g., [141 Lemma A.3]), we end up with 

^(t) = for a.e. t G (0, t) , 

whence the desired uniqueness. 

We prove estimate (|3.23|) by testing (j3.24p by "d. Using the definition (j3.ip of we find 



(3.29) 



the latter inequality ensuing from the chain rule for the convex functional j*. Hence, upon 
integrating in time we get 



(3.30) 



= l\\A\l + ^ + j/, +h + h + h. 

where the integral terms /j, i = 1,2,3 are specified by ()3.32p - (l3.34p below. Now, in view 
of (|2.H2p (see Lemma lA.ll later on), there exist Ci, C2 > not depending on e G (0,1) such 
that 

I j;{j^^.m)))>ci\\m\Wm-c2, (3.31) 

Jn 

whereas by ()3.7p and p.9p - ()3.10p we estimate the first three summands on the right-hand side 
of (|3.30p . Further, we estimate 

h= f I X^s^< /VllL6(r.)||^.||L3(r,)||^?||L2(r.) 
Jo JFc Jo 

<C f\\X\\Hiir.)0s\\mr.) m-m{^)\\v + \\mmv) 
Jo 

<C [ \\X\\Hi(r.) IWshHr.) imUnn) + C [ WXy^^r.^ WMlht.) WnLHn) (3-32) 
Jo Jo 

+ C /Vki(r.)ll^.||L3(r,)||^||Li(Q), 



Bonetti, Bonfanti, Rossi / Thermal effects in contact with adhesion 



18 



where the second inequahty follows from (j2.4p - (j2.5p . the thhd one from the Poincare inequality, 
and, thanks to Young's inequality (j2.2p . the last one holds for a suitable p > to be chosen 
later. In the same way, we estimate 



f [ div(ut)^?< / \\diY{ut)\\HmH 

Jn Jo 

<C [ \\ div{ut)\\H iW'd - m{^)\\v + WmmW) (3.33) 
Jo 

<p[ \m\h^n) + Cp' [ \\nt\\w + C [ ||ui||w!l^||Li( 
Jo Jo Jo 



(3.34) 



\v' + C fwhWv'MLHn)- 
Jo 



Now, collecting ()3.30p - (l3.3ip and ()3.32p - (l3.34|) (in which we choose p < 1/6), taking into ac- 
count (|3.26p . estimate (|3.19p for u and (|3.18p for -ds, using that fourth term on the left-hand 
side of (I3.30p is nonnegative thanks to ()2.H6p . and applying the Gronwall Lemma, we infer that 
there exists M2 > 0, independent of e and /i, such that 

e'/^\\m\\v + m\LHo,t;V) + Wrnhnn) < M2 Vt G (0,t]. (3.35) 
Secondly, we test ()3.24p by '&f Being monotone, we easily see that 

[ dt^f,{^)^t>0 a.e. in(0,r). (3.36) 
Jn 

We now integrate in time: taking into account (I3.36p . using Holder's inequality and the Sobolev 
embeddings (|2.4p - (|2.5p . we easily conclude 

4MhiO,fy) + l\Nm\\l < 1\NA\h + C mL^iT.)\\nLHT.)¥t\\L^iT.) 

+ C [ \\nL<^{rJWs\\L^rJMLHr.) 
Jo 

+ C [ \\diY{nt)\\H\\MH+ [ IMvlW: 
Jo Jo 



I.. . „o 1 



t\\v 

(3.37) 



2 

t 



< 211^011^+2^11^*1122(0,*;^) 







+ C'' / l|3<^|lHi(re) ( ll^llL4(re) + ll^sllL3{r<,) 



+ Ce" ||u,||2 + / n.n2 



Jo 

Eventually, in view of ()3.19p , p.lSp , ()3.26p , and ()3.35p , from (I3.37P and the Poincare inequality 
we deduce 

WMl^oav) + WmWv < Ml yt G (0,t], (3.38) 

from which (|3.23p follows. □ 

Due to Lemma 13.51 we are in the position of defining the solution operator associated 
with IK2^i 

72 ■■ vr2(St) X Vt ^ Wt := € HHo,t;V) : myio,t;V)nL^io,v,LHn)) < M2} . (3.39) 
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Step 3. We fix t?^, X) G Wt x 7r2(St) x 7r2(Vt) and consider (the Cauchy problem) for (fXTIll 
with data (d,ds,X). 

Lemma 3.6. Assume (l2Jn]) . (l2Jl2j) . (l2Jl6l) . (l2JT7l) . and (IgTSl) . 

Then, there exist M3, M3 > suc/i t/iai /or a// {-d, X) € Wt x 7r2(§t) x 7r2(Vt) there exists 
a unique i^s S {0,t; H^^Tc)), with 

ll^s||L2{0,t;Hi(rc))nL°°(0,t;Li(rc)) < ^3, l|(^t^s|lL2(0,t;Hi(rc)) < ^3 ' (3.40) 

sitc/i i/iai complies with the initial condition (|3.16|) . and 

ej dt^sv+ I dt^^{ds)v- I dt\{X)v+ I V^.Vv (3.41) 
+ e /" Vdt^s'^v= [ k{X)(d-ds)v yveH^{Tc) a.e. m(0,t). (3.42) 

Proof. Thanks to |17l Thms. 1,4] there exists a unique solution -ds € H^{0,t;H^{Tc)) to 
the Cauchy problem for ()3.4ip . Hence, we shall just prove ()3.40p . referring to notation (13.250 
for the term k{X). We proceed as for (|3.23|) : hence, we test (|3.4ip by i!)s and integrate in 
time. Developing the very same calculations as throughout (I3.29p - ()3.31|) . exploiting (I2.H2P (via 
Lemma EH), and recalhng I^B) and (l3:9]l - (|3l^ . we find 



e ,, „ , 1,9 ^ I, „ . .X 11 ; ,,„ „ ,|2 



+ Ci\\Ut)\\LHr.)+ / l|Vi?.|| 



2" 

< C + /4 + /5 + ^6 







where, also in view ()2.H7p . 

ft r ft 



h= l^t||A'(X)||T?,| <C / ||Xt||i2(r,) (||X||i4(r,) + l) ||^,||i4(r,), (3.43) 
Jo JVc JO 

and, in view of (|3.26p . 

miHr.) PllL4(r.) IWsWlhf.) < l|3<:||L-(o,t;L4(r.)) / Plk^cr.) WMl^f.) , (3.44) 

Jo 



h < f 

Jo 



h= j / P^\\L^(r^)\\^s\\L3(rc)\\'^s\\L^{rc) 



° -^^^ ^ (3.45) 



< l|3<^llL°°(0,t;L6(r,)) / \\'^s\\l3(V,)\Ws\\l^{V,)- 

Jo 

Taking into account (I3.18|) . (I3.22|) . and ()3.39p . and estimating the term ||'!?s||L2(r^) by ||i9s||Li(r;,) 
and II Vi?s||^2(-p ) in the same as in (|3.32p - (|3.33p . we finally apply the Gronwall Lemma to 
conclude 

£"^''^Ps||L°°(0,t;Hi(rc)) + ll^s||L2{0,t;Hi(rc))nL°°(0,t;Li(rc)) - ^3- (3.46) 

Then, we test (13.41 p by dt'&s and integrate in time. Thanks to ()3.46p and arguing in the 
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same way as in (j3.43P " (|3.44p . we find 

e fWsfH^(^r.) + l\\Mt)\\mir.)+ f dt^,{^s)dt^s 
1 /■* _ _ 

<^¥"s\\mir.) + \\M\L^(V.)\\nLHT.)\\dt^s\\L^iV.) 

+ f \\nLHV.)¥\\L^{V.)Ws\\L^(T.) 

Jo 

+ [\\nLHr.)mLHr.)Ws\\LHr.)- 
Jo 

Again taking into account (I3.18|) , (I3.22|) , and (I3.39|] , as well as the fact that the third summand 
on the left-hand side is nonnegative, we readily deduce the second part of (I3.40p . □ 

Lemma 13.61 enables us to define a solution operator associated with (|3.4ip 

Tg : WtX7r2(§t)x^2(Vt) ^ := {^s G 1°^ {0,t; H^T,)) : \m\LHo,v,mir.))nLoo^o,v,LHr.)) < Ms} • 
3.3 Local existence for Problem (P^) 

Proposition 3.7. Assume (12111]) - (|2.H10I1 . (l236l) - (ICTH) and (IHTp - ffXTOD . Then, there exists 
T € (0, T], possibly depending on fi > 0, such that for every e > Problem (Pe) admits a 
solution (i?, u, X, ^, 77) on the interval (0, T). 

Proof. In view of the above Lemmata l3.4l [33tl3.6l we are able to define an operator T whose 
fixed points are solutions of Problem (Pg ). 



Definition of T. In the end, we define 

T:St^WtxytX7r2(V0 

by setting for every i?s,X) S St 



T(t?,t9„X) := where < 



^s = 'Js{M^s, 7i{d,ds,x)), ds, 7r2{7i{d,ds,m, (3.47) 



Thus, in order to prove Proposition 13.71 it is sufficient to show that there exists T G (0,T] such 
that for every e > 

T maps into itself, (3.48) 
T : — > is compact and continuous 

(3 49) 

w.r.t. the topology of L^{0, f; H) x L^I^{Q, T; L^Tc)) x L^°{0, f ; L^{Tc)). 



Ad dOH]) . We fix {'d,'ds,X) G St and let (t?,i?s,X) = T(i?, i?,, X). First of all, the three- 
dimensional version of the Gagliardo-Nirenberg inequality (cf. |23^ p. 125]) yields 

||^k<C||^||f ||^||^(^(^), (3.50) 
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so that, using ()3.23p . we get 

ll^llL^(o,t;^.) < m%%,,,v) \\n%%,,,L^n)) < . (3.51) 

The Gagliardo-Nirenberg inequahty in 2D gives 

||^.||L3(r.)<C||^.|lH/iV)ll^.lli((rc)' 

whence 

ll^.llL5/.(o,t;L3(r.)) < m\%%,,HHr.)) ll^lli^(o,t;Li(r.)) < ''^''M,. (3.52) 
Finally, by (13.190 we have 

ll>^llLio(o,t);L6(r.) < Ct^/iOMi , (3.53) 

the constant C depending on the Sobolev embedding ()2.5p . Clearly, there exists T (which does 

not depend on e) such that (i?, t?s,X) belongs to 8^, hence the operator T maps into itself. 



Ad (I339|). Exploiting (IXT9|) . (I3:23D . and (^M), the Sobolev embeddings (l23D-(l23]) and 
Thm. 4, Cor. 5], one sees immediately that the operator T : 8^ — > is compact. We shall prove 
that T is continuous in three steps, basically checking that the operators 7i, i = 1,2,3 defined 
in Section [32] are continuous w.r.t. to suitable topologies. 

We fix a sequence {{'&m''^s,n,^n)} C §f such that there exists (t?oo) i^s.oo ^cx>) S with 

as oo, (3.54) 
^s,n^^s,oo in L^^\0,f-L^{Tc)) as n ^ oo , (3.55) 
^n^^oo in L^°(0,f ;L^(rc)) as n ^ oo , (3.56) 

we let (u„,X„) := 7i{'dn,'&s,n,^n), and denote by {t7,j} and ^„ the associated sequences of 
selections of the graph a and /?, respectively, such that (j2.37p and (j2.39p hold for all n G N. Due 
to (I3.19p . there exists a subsequence (which we do not relabel) and a quadruple 

(uoo,r7oo,^oo,eoo) e H\0,f;W) x L\0,f;H~^/\T,)) 

X (L2(0,f ; F2(re)) n L°°(0, f ; //^(re)) n i/HO,T; L\T,))) 
xL\0,f;L\r,)) 

such that the following convergences hold as n oo: 

u„-Uoo in H\0,f;W), u„ ^ u^o in C'^{[0,f];{H^-'{n)f) for all s > 0, ^ ^ 

n - , (3.57) 
u„ ^ Uoo in C°([0,r]; (^^(re))^), for ah 1 < p < 4, 

Vn-Voo in L2(o,f;(F-i/2(r,))3), (3.58) 

X„-*Xoo in /7i(0,f;L2(r,))nL°°(0,f;Fi(r,))nL2(0,f;F2(r,)), 

^ ^ (3.59) 
X„^Xoo in C°([0,r];//i-^(rc))nL2(0,r;F2-^(re)) foralls>0, 

e„^^oc mL\0,f;L\rc)). (3.60) 

Thanks to these convergences, and recalling ()2.H7p . it is possible to pass to the limit in the 
PDE system (|3.2Up as n /' oo (see the proof of [U Prop. 4.7] for further details). Thus, we 
prove that the quadruple (uqo, ■/7oo; -^oo, Coo) fulfils ()3.20p . with data {'&oo,^s,oo,^oo), on (0,T). 
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Furthermore, the pair {uoo,^oo) clearly complies with the initial conditions ()2.3ip - (|2.30p . Note 
that the identifications 

Coo e /3(Xoo), Voo e "(Uoo) 

are proved by semicontinuity arguments (see O Lemma 1.3, p 42]). In the end, recalling the 
definition of the operator Ti we conclude that 

(Uoo,Xoo) = Ti(l?oo,l9s,oo,Xoo) • (3.61) 

In fact, by uniqueness of the limit we readily deduce that convergences (j3.57p - ()3.60p hold for 
the whole sequences. In particular, we have checked that (|3.54P " (|3.56p imply 

'Ji{dn,ds,n,Xn) ^ Ti (i9oo, ??s,oo, ^oo) in the scnse of (I33ZD, ^M- (3.62) 

We now consider the sequence '&n := 72{'ds,n,'^n,^n) = '^2{'&s,n,7i{'dn,'ds,n,^n)) ^or all 
n e N. Thanks to (13:23]) . is bounded in H'^{0,f,V), hence there exists t?oo S H^{0,f,V) 
such that (up to a subsequence) 

i?„-i?oo in H\0,f-V), ^n^^oo in C°{[0,f];H^~%n)) Vs > 0, 
^oo in C°([0, f];LP{rc)) for all 1 < p < 4, 

as n ^ oo. In particular, "i^co complies with ()3.15p . Moreover, using that is Lipschitz 
continuous, taking into account the strong convergence for i?„ specified by the second of (|3.63p 
and recalling that, by maximal monotonicity, the graph of is strongly- weakly closed, it is 
not difficult to conclude that 



^^{^n)^*^^.{^oo) in L-(0, f; V) n H^O, f ; H), 
^^{dn) ^ ^^,{^oo) in C0([0, f]-H). 



Furthermore, it follows from (|3.55p . (|3.59p . (j3.63p . and the Lipschitz continuity of k that, among 
others, the following convergences hold n /' oo: 

KXn)^n ^ A;(Xoo)i9oo in L°°(0, f ; L2(rc)) for all 1 < p < 4, and ^ ^ 

^ - o (3.65) 
k{Xn)'&s,n ^ A;(Xoo)i?s,oo in l2(0, T; L^iT,)) . 

Combining ()3.57p with ()3.63p - (l3.65p . we pass to the limit as n oo in (I3.24p with data 
{'ds,n-, Un, Xn), and we deduce that ^oo solves (the Cauchy problem) for equation (j3.24p . with the 
triple (t?<i,oo,Uoo,Xoo), on (0, T). Thus, 

f^oo — '^2 {,'&s,ooi Uqo ; '^oo) — '^2 (^s,oo j '^1 i'&oo ; ^s,oo i -^oo)) ; 

the second equality ensuing from (I3.6ip . Again, since the limit t?oo does not depend on the 
subsequence in (j3.63p . it turns out the convergences specified therein hold along the whole 
sequence {i^n}- In conclusion, 

'^2{'&s,n,7i{dn,'&s,n,Xn)) ^72{ds,oo,7i{doo,'&s,oo,Xoo)) in the sense of (|3.63p . (3.66) 

Finally, we let 

^s,n := 73{^n,^s,n,Xn) = 73{M'^s,n, Ti (^„ X„)) , TTs (Ti (l?„ X„)) ) 
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for every n G N. By (j3.40p . {'i?s,n} is bounded in H^{0,T; H^{Tc)). Thus, there exists a (not 
relabeled) subsequence and ??s,oo ^ {0,T; (Tc)) such that 

^s,n ^^s,oo mH\0,f;H\r,)), ^ 

n - 1 (3-67) 

^s,n ^ ^s,oo in C°([0,r];/7i-^(r,)) for all s > 0, 

^^.{^s,n)^*^f.{^s,oo) mL^{0,f-H\r,))nH\0,f;L\r,)), 

if^(i?,,„) ^ if^(i?,,oo) inC0([0,f];L2(r,)), 

where the latter convergences can be proved arguing in the very same way as for ()3.64p . In par- 
ticular, i)s,oD fulfils initial condition (j3.16p . and, like in the proof of Lemma 13.61 combining con- 
vergences (|3.67p - (j3.68p we conclude that '&s,oo fulfils (j3.4ip on (0,T), with data {'&oo,'&s,oo,^oo)- 
Now, (I3.59P and the growth properties of A (see (I2.H7P ) clearly yield that 

X'{Xn)dtXn^ X'iXM^oo mL\0,f;H\T,y). (3.69) 

Then, exploiting (|3.59p and (j3.67p we easily check that 

KXn)i}s,n ^ A:(Xoo)i?s,oo in L°°(0, f ; L^T,)). (3.70) 

Collecting (|3.67p - (j3.7Up and also taking into account (|3.65p . we pass to the limit in (|3.4ip (with 
data {'&n,'&s,n,^n)) as n y oo. Hence, 

,005 '^2('J'l(^00) ''^S,005 -^oo)))) 

and we again deduce that convergences (I3.67p - ()3.68p hold for the whole sequences {i}s,n}- Even- 
tually, we have that 

T3(T2(T?,,n, 7i{dn,ds,n,Xn)), ^9s,n, VTz (Ti (l?„ X„) )) ^ 

^ ^ ^ ^ ^ ^ ^ ^ {o.iL) 

'J'3('J'2(i9s,oo, Ti{'&oo,'&s,oo,Xoo)), '&s,oo, 7r2('Ti (t?oo, i9s,oo, ^00))) in the sense of (|3.67p . 
Clearly, (j3.62p , (j3.66p and (|3.7ip show that T is continuous in the sense of (j3.49p . n 

3.4 Global existence for Problem (P^) 

In order to show that the local solution to Problem (Pe ) actually extends to the whole time inter- 
val (0, T), we shall obtain some global in time estimates on the solution components ("!?, u, X) 
and then use a fairly standard argument to conclude that, for every e > and /i > 0, the local 
solution found in Proposition 13.71 extends to the (unique, by the calculations in Section 13. 5p 
global solution of Problem (Pg ). 

Lemma 3.8 (Global estimates). Assume (|2.IIip - ()2.II10p and let (■!?'', -iJg, uq, Xq) be a quadruple 
of initial data complying with conditions ()3.7p - 1)3. lOp and ()2.16p - ()2.17p . Then, for every fi > 
there exists a constant > 0, depending on the problem data, on Mq ( cf. (j3.9p ) and Mq 
(cf. ()3.10p ), but neither on t G (0, T] nor on £ > 0, such that for every solution (i?, t?^, u, X, rj) 
to Problem (Pe ) on the interval (0, t) there holds 

^^^^\\'&\\L°°{0,t;V) + ll'*9||L2{0,t;V)nL°°{0,t;Li{Q)) + II J/^ (-^M (^)) (0,t;Ll (f^)) - ^/^ ' (3.72) 
^^^'^\Ws\\l°° {0,t; H^Tc)) + ll^^s||L2(0,t;/fl{rc))nL°°(0,t;Li{rc)) 

+ \\j;{^M)\\L^mLHr.))<C^, 
ll^llHi(o,t;L2(r^))nL°°{o,t;Hi(rc)) ^Cf,, (3.74) 
l|u||//i{o,t;W) < C/j,. (3.75) 



(3.73) 
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Proof. We test ([2361) by uf. owing to and to (f2ll3l) . via the chain rule for a and 

the Holder inequality we obtain 



^ / l|ut|lw + ^l|u(t)||w+ / i?div(ut) + a(u(t))+ / / Xu ut 

^ JO ^ JQ Jo JTc 

<C||uo||^ + a(uo) + ^^ ||F||2.,. 



(3.76) 



Next, we multiply (j2.38p by Xt and integrate in time. Again applying the chain rule to the 
functional and using (|2.H5p . with easy calculations we find 



1 r 1 



<C + ^llVXo||i2(r^)+ / mo)- [ [ \'{X)Xtds + l7 + h. 



(3.77) 



where 



h = Llj'^j^ |X|2<2L2r||Xok2(r^) + 2L2r^*(^^'||Xt||2,(^^^^ ds (3.78) 
while, with an easy integration by parts, we have 

h = -\l I M^?=! I ^^tvi-l [ + toluol'. (3.79) 

Furthermore, being /? convex, we have 

Mm) > -Ci,f3Mt)\\mr.) - C2,p 

>-v\\m\\h^r.)-C^ (3.80) 

> -2??r^ l|Xdli2(r,) - 2??||Xo||i2(r^) - C, 

for some suitable r/ > to be specified later. Finally, we test (j3.13p by i}, (|3.14p by i?^, integrate 
in time and develop in both cases the same computations as throughout (I3.29p - (l3.3ip . We add 
the resulting inequalities with (j3.76p and (|3.77p : also taking into account (j3.79p . some terms 
cancel out. Furthermore, choosing ij < 1/8T in ()3.80p and applying the Gronwall Lemma to 
deal with the integral term on the right-hand side of (j3.78p . we arrive at 

|2 , ||„Q||2 , |i„a^+MI . , ^ILQ ^+M|2 , ||„Q ||2 



{o,t;V) + II^WIIlmh) + e|l^s(t)|lHi(r,) + ll^s|| 
+ ll^s(t)||Li{r.) + ^ km^-^s)' + \\MUo,t;W) + Mt)\\w 

+ [ mHt)\^ + «(u(t)) + \\M\h^o,t;LHr.)) + l|VX(t)||i2(r,) 



(3.81) 



<c(i + n\h\\i,+ I iiF"2 



w 



for every t E [0, T]. Noting that the tenth and eleventh summands on the right-hand side are 
nonnegative due to (|2.H3p and (|2.H4p . we conclude. Ultimately, we also find that 

Wf.i^.mh^mLHn)) + \\j;{-^,{^s))\\L^io,v,LHr.)) <C, (3.82) 
and (l372]) - (f3:75]) ensue. □ 
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Remark 3.9. As it is clear from the above proof, constant in (j3.72p - (j3.75p depends on 



only through Mq (j3.1(jp . i.e. the bound for and £^^'^\\Vs,eiJ.\\H^{rc)- 

3.5 Uniqueness for Problem (P^). 

We prove here the uniqueness statement in Theorem l3.ll Let us consider two families of solutions 
"iJs^i, Uj, Xj, .^i, ?7j), i = 1,2, to Problem (Pe). Hereafter, we shall refer to the notation 

^='di-'d2, -ds = '&s,l - ^s,2, U = Ui-U2, X = Xi-X2 ^=^1-^2, V = - V2- 



We will derive suitable contracting estimates on the solutions. First, we subtract (|3.13p written 
for (i?2 5 i9s,2, U2, X2) from (I3.13P written for (■!?!, -i^g^i, ui, Xi) and integrate in time. We get 



e / d{t)v+ / {^^{di{t))-^^{'&2{t)))v + e / W^{t)Vv+ / {l*\/^)(t)Vv (3.83) 
Jn Jo. Jn Jn 

= divS(t) V + (l * [k{X2){^2 - ^s,2) - h{Xi){^i - l9,,i)]) (t) V 

for all V ^ V . Letting v = "d, integrating in time and exploiting the monotonicity of .if^, we 
obtain 



Jo ^ Jo Jn 



(3.84) 



In order to estimate the three latter summands, we proceed as follows. Using Holder's and 
Young's inequality (cf. also ()2.ip ). and well-known Sobolev embeddings, we find 



< 



Jr. 



L2 



l*[ik{X2)-k{X^))i^l-^s,l)]y 



(r,) ||1* [(A;(X2)-A;(Xi))(^i-^?,,i)]| 



lL2(r,) 



< C 







IlHt.) \\{K^2) - k{Xi)){^i - ^s,i)\\mo,s;L^r.)) ds 



(3.85) 



< C / ms)\\L2^r.) \\kiX2) - A;(Xi)||i2(o,,;L4(r,))||i9i - ^?s,i||L-(o,s;L4(r.)) ds 
Jo 



< 5 



for a suitable positive 6 to be chosen later. In particular, the positive constant cs^i also depends 
on the a priori estimate (|3.72p on — '&s,i\\L°°{o,T;L*{rc))- Arguing similarly, also taking into 
account (j3.74p in order to estimate X2, we have 



'10 



Jr. 



1* [A;(X2)i?])i? 



< 







lL2{r,) 111 * [k{^2)'&] llL2(r,) 



<C ||i?(s)||L2{r,) 11^2 + l||L-(o,s;L4(r,))||i9||L2(o,,;L4(r^)) ds 
Jo 

<S I ll^?l|^ + C5,2 / ll^?lli2(o,.;y)ds. 



(3.86) 
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and 



'11 



< C 



< 6 



Jr 



1 * [k{X2)^s\)^ 



< 



L2(r,) 111 * [^(^2)i9s]||L2(r,) 



(■5)||L2(r^) \\X2 + M\L°°(0,s;L'^(rc))\Ws\\L2{0,s;L*(rc)) 

ft 



(3.87) 



V 



+ CS,3 



^slli2(o,s;Hi(rc)) ^•^ ■ 



In a similar way, we consider ()3.14p written for [1^2, "& 8,2,^2) and for (i?!, Xi), we take 
the difference, test the resulting equation by t?^ and we integrate in time. We obtain 



(3.88) 



(re) 



where the terms I12 and /13 are estimated in the following way. Recalling (j2.1ip . we find 



(12 



Jr. 



A(Xi)-A(X2)U 



<C / ||X||i2(r,)(||Xi||i4(r,) + ||X2||L4(r,) + l ||^9s||L4(r.) 



(3.89) 



< 6 



+ CSA 



\x\ 



ds. 



where the constant 5 > is the same as in (I3.85p -( r3.87p . Note that the positive constant 05^4 also 
depends on the estimate for ||Xi ||j;^oo(o,r;L4{rc)) ll-^2||L°°{o,T;L4(rc)) given by (j3.74p . Moreover, 
arguing as in the derivation of (|3.85p - (|3.87p . we have 



'13 



Jr. 



1 * [k{Xi)i^i - - k{X2){^2 - ^s,2)\)^ 



X 



L2{0,s;//i(rc)) 



(3.90) 



+ C5,6 / l|l?llL2(o,,;y) +C5,7 / ||l9s|lL2(o,,;/^i(r,))ds. 



Now, we subtract (j2.36p . written for (i?2, U2, X2, ^72)) from (|2.36p . written for (i?i, ui, Xi, rjj^), 
we test the resulting relation by u and integrate on (0, i). Recalling (|2.7p - (|2.8p . and the mono- 
tonicity of a (cf. (12.H3P ). we end up with 



a, 



w 



+ Ca||u| 



< 



L2{0,i;W') 



X2 (u) 



+ 



Jr. 



diviJu 



JTc 

t 



Jrc 



Xuiu< ||X||i2(r^)||ui||i4(r^)||u||i4(r^) (3.91) 



< C I IIXll?^,^ ^ + C / IISII2 



iL2(r,) 



W ' 



where the last two inequalities follow from Holder's and Young's inequalities, from Sobolev's 
embeddings and from the fact that X2 > a.e. on (0,T) X Tc, due to (ITHi]) . In particular, the 
constant C in (j3.9ip depends on ||ui ||2,oo(o.T;L4(rc)) through estimate (|3.75p . 



Bonetti, Bonfanti, Rossi / Thermal effects in contact with adhesion 



27 



On the other hand, let us consider the difference of (|2.38p written for ui, Xi, ^i) and 
(imSI) for ('i?s,2,U2,X2,C2), multiply it by X and integrate on (0,t) x r^. Taking (I2ll4j) . (l2TH5l) . 
and (I2.H7|) into account, we get 



t 



X'{Xi)-X'{X2))^s,2X-l I I (ui + u2)-nx 
< L^l^ \\X\\hir.)+cl^ (l|Xi||L4(r,) + l)\Ws\\LHr.) MlHt.) (3-92) 

+ C \Ws,2\\L^{rc) ll^llL'i(r^) ||^||L2(r^) + C ||ui + U2||/^4(r^) ||u||i4(r^) ||^||L2(r^) 
Jo Jo 

<5 ¥s\\]mr.) + 5 \\xfmiT.) + cl \\n\\w + cs,8j^ mh^r.), 

where the constant depends on estimates (I3.73p - ()3.75p for the quantities ||Xi ||j^oo(o,T;L4(rc))5 
||'(?s,2||L°°(o,T;L4(rc))'^;nd onjlui + U2\\L^{0,T-,lA(rc))- 

Finally, we add ()3.84p , (I3.88|) , (I3.91|) , and ()3.92p . Noting that two terms cancel out and taking 
into account (|3.85p - (|3.87p and (|3.89p - (|3.9Up (in which we choose, e.g., 6 = min{e/6, 1/2}), we 
find ^ ^ 

^ ~^ ^ Jo ''^*''^^(^':) ll-^ll/fi(re) + Il-^(*)lli2(rc) + ^^'ll"(*)llw 

~^ Jo "^"^^(O'^^^^rc)) ^'^ + ll^llL2(rc) + ^ l|u|lw 

Thus, by Gronwall's Lemma, we conclude that '&i = '&2, '&s,i = '&s,2i ui = U2, and Xi = X2. 
A comparison in (j2.36p and (j2.38p also yields iji = 772 and ^1 = (,2, so that the uniqueness 
statement in Theorem |3JJ follows. □ 



4 Proof of Theorem [T] 

As we mentioned in Remark 13. 2| we shall pass to the limit in Problem (Pe) first as e \ and 
/X > is fixed (cf. with Proposition 14. 4p . and then as /x \ (see Section Iir2]) . The next result 
(whose proof is postponed to the Appendix) concerns the construction of sequences of initial 
data C V and {1?°^^} C H^{Tc) for Problem (P^) complying with (l371l- (l3lnD . and such 

that the sequence of solutions to Problem (Pe), supplemented with the data (-i?^^, -iJ^^g^, uq, Xq), 
converges to a solution of Problem (P) in the two consecutive limit procedures e \ and \ 0. 
As, in our construction, the data in fact depend only on the parameter > 0, we shall denote 
them as and 1!}^^^ for simplicity. 

Lemma 4.1. Assume that the initial data wq and zq respectively comply with ()2.14p and (|2.15p . 

Then, 

1. there exists a sequence {uP^}^ C V fulfilling for every fi > 

\\wf,\\H < \\wo\\h , (4.1) 

/ f < / f (wo) , (4.2) 
Jn Jn 
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and such that 

—>■ Wo in H as /i \ 0. (4-3) 

Furthermore, let us set 

^l:=^-\wl) for all ^i>^. (4.4) 

There exists a constant C^o > 0, depending on wq and on fi > but independent of e > 0, 
such that for all e > 

e'/'\Wl\\v<C^:^,. (4.5) 
2. There exists a sequence {z^}^ C //^(Fc) fulfilling for every fi > 

lk°llL2(re) < lko||L2{r,) , (4.6) 

/ f {z',) < [ f izo) , (4.7) 

and such that 

zl zo in L^{Tc) as fi\0. (4.8) 

Furthermore, setting 

:=^m"'(^°) foraUf,>0, (4.9) 

there exists a constant Czq > 0, depending on zq and on n > but independent of e > 0, 
such that for all e > 

e'^'\\€jmir.)<C^,- (4-10) 

Remark 4.2. Note that our construction of the sequences and only depends on the 

data wq and zq. This is, ultimately, the main reason why we have chosen to approximate the 
operator £ by the regularization .i?^, instead of the usual Yosida regularization Indeed, if we 
had used the latter approximation of i, starting from the datum wq we should have constructed 
approximate data "d^ G V satisfying the corresponding bound 

WUf^nWlmLHU) < C for all ^ > 0, (4.11) 

(where (j^)* is the conjugate of the Yosida approximation of j), cf. with the proof of 
Lemma 13.51 (clearly, the same considerations hold for i^^^^J. To deduce ()4.1ip from the con- 
dition fiwo) G L^(Q), one should virtually choose i)^ in such a way that ^/^(i?^) = wq (leaving 
aside the condition € y). However, it is not clear to us how to carry out this construction, 
since is not invertible. Instead, .5?^ can be inverted, and the calculations we shall provide in 
the proof of Lemma l4. 1 1 show that the sequence defined by (14. 4p complies with (I3.9p - ()3.10p . 

Notation 4.3. We shall denote by 

{{^e^l,'&s,e^l,'^e^l,^e^l,^e^l,1^etJ.)} the sequence of solutions to Problem (P^) 
with initial data {(t?^, t?^^^, uq, Xq)}. 

Further, for simplicity we shall use the notation 

SO that, in view of (|4.4p and of (|4.9p respectively, we have 

«;,^(0) = u;°, 2,^(0) = z° for all e > 0. (4.12) 
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4.1 Passage to the limit in (P^) as s \ 

Proposition 4.4. Let /_i > 6e fixed. Under the assumptions of TheoremUl there exists a (not 
relabeled) subsequence and (t?^, w^, -i^^^^, z^, u^, X^, 77^) such that the following convergences 
hold as e \ 

^e^.^^f. inL\0,T;V), eJl{^ei.) ^ in L°^{0,T;V'), (4.13) 

^s,ef.^^s,^. inL\Q,T-H\T,)), e-RrMs,e^.) ^ ^ in L^{0,T; H^T,)'), (4.14) 

We^^Wf, inL\0,T-V), (4.15) 

eJl{^s,.)+We^^Wf, in H\0,T;V') , (4.16) 

Ze^^z^ tnL\0,T;H\r,)), (4.17) 

eJl^M,e^.) + Ze^. ^ z^ tn H\0,T; H^F,)') , (4.18) 

Xef,^*X^ in H\0, T; L^T,)) n L°°(0, T; H^F,)) n L^{0, T; H^F,)), 

X.M^^M inC\[0,T];H^-^{F,))nL\0,T;H^~\F,)) for all 5 > 0, 

C,^-^^ mL^{0,T;L^{F,)), (4.20) 
2, 



(4.19) 



Vs^^Vf. tnL\0,T;W'), (4.21) 
u,^^u^ inH\0,T;W), 

Usi, u^, m C^{[0, T];H^~^{nf) for all S >0 . 

Moreover, and z^ have the further regularity 

G L°°(0, T; H), j; (w^) G L~(0, T; ^^(f])) , 
G L~ (0, T; L2 (r,) ) , j;{z^)eL^{0,T;L\F,)), 



(4.22) 



(4.23) 



(4.24) 



anc? satisfy 

f u;^(x,t) = Jf^(i?^(x,t)) /or a. e.(x,t) G X (0,r), 

\ Zf,{x,t)=^f,{^sA^^t)) for a.e.{x,t)eFcX{0,T). 

Further, the functions {■d^,w^,'ds,^i, z^,Uf^,X^,^^,r]^) fulfil equations (|2.36p - (j2.40p and (|3.13p - 
p.l4p with e = 0, and the quadruple {w^^z^^u^^X^) complies with the initial conditions 

w^{Qi)=wf, a.e. inn, z^(0) = z° a.e. inFc, ^^^^^ 
u^(0) = uq a.e. in J7, X^(0) = Xq a.e. in Fc. 

Notation 4.5. Hereafter, we shall call (P^) the boundary value problem given by ()2.36p -( r2.40p 
and (j3.13p - (|3.14p with e = 0, supplemented with relations (j4.24p and the initial conditions (j4.25p . 



Proof. Let us point out that, thanks to (|4.ip - (|4.5p and (|4.6p - (|4.7p . all estimates in Lemma [ 
hold for the sequence {i'^e^J.,We^l, '&s,eij., ^e/^, u^^, ^e/^, CeiJ.,'nefi)}e with a constant only depend- 
ing on the problem data, on the initial data {wq, zq,uq,Xo), and on /i > 0, but independent of 
e > 0. In order to pass to the limit in (Pe ) as e \ 0, we need some further estimates in addition 
to ()3.72p - (j3.75p . Using the latter bounds and arguing by comparison in (j3.13p and in (j3.14p . we 
conclude that there exists C > such that for all e, /i > 

\\e3l{dtT9s^) + dtWsjL2(o,T;V') + Ik3?r,(<9ti?s,e/.) + dtZejL^o,T;m{r,y) <C . (4.26) 
Similarly, a comparison in (j2.36p leads to 

N£mIIl2(o,T;W') < C for all e, /i > 0. (4.27) 
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Finally, we test (j2.38p by E Pi^efi) and get by standard arguments 

Ue4L^0,T;L^{r,)) + IIl2{0,T;//2 (Fe)) < ^ for all E, fl > 0. (4.28) 

Moreover, in view of the Lipschitz continuity of .if^ (cf. with ()3.6p ). joint with estimates (I3.72p 
and (|3.73p for and i?^, we conclude that there exists a constant > 0, depending on the 
problem data and on /i > but independent of e > 0, such that 

\\'Weti\\mO,T;V) + lkeMllL2{0,T;Hi(rc)) < C'm £ > 0. (4.29) 

Combining estimates (|3.72p - (j3.75p and (|4.27p - (|4.29p with the Ascoli-Arzela theorem, the well- 
known \28\ Them. 4, Cor. 5], and standard weak compactness results, we find that there exists 
an eight-uple {'&fj_,w^,'ds,fi, z^,u^,X^,ri^,(,^) such that, along a suitable (not relabeled) subse- 
quence, convergences (Iil3]) - (lil5]) . (filTp . and (Iil9l) - (|i:22l) hold. Clearly, (filH]) follows from 
(j4.26p and the second of (j4.13p . In the same way, we obtain (j4.18p . Therefore, the further regu- 
larity ()4.23p for and ensues from the continuous embeddings L^{0, T; V) fl H^{0, T; V') C 
C^{[0,T];H) and L^{0,T; H^Tc)) n H\0,T; H^iTc)') C C^{[0,T]; (Tc)). In order to prove 
the second of ()4.23p . we exploit a Lebesgue point argument. Indeed, by ()4.15p and the lower- 
semicontinuity of the integral functional induced by (the convex function) j* with respect to 
the weak convergence in L'^{0,T;V), we find that for all to ^ (0)^) and r > such that 
{to-r,to + r) C (0,T) 

/ / j;(.w^) < Hminf / / j;{ws^) < 2r sup \\j;{w,^)\\L^io,T;LHn)) < 2rC , (4.30) 

the latter inequality due to ()3.72p . We now divide the above relation by r and let r | 0. Using 
that the Lebesgue point property holds at almost every tQ G (0, T), we obtain the estimate 

lb;K)llL-(o,T;LMn)) <C^ forall/i>0, (4.31) 

C being the same constant as in (I3.72p . The analogous bound for j* (z^) is proved in the same 
way. 

Furthermore, we remark that ()2.H6p . ()2.II7p . (I4.13p ^ ()4.19p . trace theorems, and Sobolev 
embeddings yield that, as e \ 0, 

k{Xe^){^,^ - 1?,,,^) - k{X^){^^ - ^s,^) in L2(0, T; L^{T,)), 

(1 * kiXe^){^e^. - ^s,ef.)) ^ (l * k{X^W^ - ^,,^)) in C^{[0,T]; H-'/\T,)), '^^^ 
as well as 

A(X,^)-A(X^) mH\0,T;H\r,y). 

In order to pass to the limit in Problem (Pe ) as e \ 0, we use the above convergences and 
proceed as in Section [331 In particular, for i^M^-l^iM) we exploit (fil3]) . (fil9]l - (|i:22]) and 
argue in the same way as in [H Prop. 4.7], to which we refer the reader. Further, relying on 
([il^ - ifil^ and the above we also pass to the limit as e \ both in ([TTH]) and 

in (I3.14|) . We thus conclude that (i?^, tt;^, i?^ ^, z^, u^, X^, ?7^, ,^^) satisfies the PDE system given 
by ^M-^M ^B-^M with e = 0. 

Finally, it remains to show (I4.24p . We shall just prove the relation for u;^, the argument for 
being completely analogous. By maximal monotonicity of (cf. Lemma 1.3, p. 42]), it 
is sufficient to show that 

limsup [ [ We^i^et,< [ [ w^"&^. (4.34) 
£\o JO Jn Jo Jn 
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which we prove integrating in time p.l3p . testing it by i^g^, and again integrating on (0,T). 
Hence, we get 

hmsup / / Wsfj.'dsiJ. 

e\0 Jo JQ 



<hmsup / / + -div(uo))i9eM + limsup / [ eV^y^e^^ - elimmf [ 

e\0 Jo Jn e\0 Jo Jo. "^^0 Jo 

+ Yimsnv f [ div(u,^)T?,^ - hminf / |(1 * V^9,^) (r)p 
£\o Jo Jn Jn 

-hminf / / (l * /c(X£^)(t?£^ - '(Js^e^))??^^ + hmsup / (l*/i,i?£^) 
e\o Jo JFc e\o Jo 



Jo Jn 

where the second passage follows from (14. 5p and the first of (14.130 . from the second of ()4.13p . 
from (j2.32p and the following relations: 

lim / {l*h,i^,f,)= [ {l*h,^^), lim / / div(u,^)i?,^ = / / div(u^)i?^ 
£\oJo Jo £\ojo Jn Jo Jn 

(due to (fil3]) and ([02]) . yielding in particular that div(u£^) div(u^) in L°°(0, T; ^^/^(O)) 
as e \ 0) , 

hminf/ |(1*V^,^)(T)|2> [ \{l*V^,){T)f 
=\o Jn Jn 

(by the lower semicontinuity of the norm, since 1 * V?9e^^*l * Vt?^ in L°°{0,T; H) as e \ 0), 
and, finally, 

lim r / {l*k{Xe^,){'&,^-^s,e^^)^e^. = r I (l * A:(X^)(i9^ - 
<^\^^Jo Jtc Jo Jtc 

thanks to the second of and KT^ (which gives that ^ i^f, in L^{0,T; H^/^(Tc)). 

In the end, we show (I4.25p . The initial conditions for and ensue from ()4.3p and 
convergences (|4.19p and (j4.22p . On the other hand, thanks to (j4.5p there holds 

eJli'd^) ^0 as e \ in y, 

hence by ([iT2]l and ([iT6]) . 

= lim u;,^(0) = lim (e3i(i?0) + We^iO)) = w^{0) , 

e\0 £\0 

where all of the above limits are meant with respect, e.g., to the -ff^(r2)'-topology (in fact, to 
the topology of any space Z such that V' C Z with compact embedding). With a completely 
analogous argument, we prove the initial condition for as well. □ 
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4.2 Conclusion of the proof of Theorem [T], 

We are now going to show that the sequence {(ti)^, -i?^, z^, t?^^^, u^, X^, 77^)}^ of solutions 
to Problem (P'^) (cf. with Notation I4.5P obtained in Proposition 14.41 admits a subsequence 
converging as /i \ to a solution of Problem (P). To this aim, we point out that there exists 
a constant C > 0, independent of > 0, such that 



II^At||L2(0,T;y)nL°°(0,T;Li(f^)) + IIJa' (-^M (^A*) ) II (0,T;L1 (f^)) - ^' 
\Ws,^l\\L■2{0,T■,H^^c))nL^{0,T■,L^^c)) + \\jl.i-^f^i'^s,fi))\\L°°{0,T;L^rc)) ^ C", 
II^At|li^M0,T;L2(rc))nL°°{0,T;_H'i{rc)) ^ ^' 
^ ll%lliTi(O.T;W) < C. 



(4.35) 



This can be proved by testing (IXTH|) (with e = 0) by 1?^, ([XTi]) (with e = 0) by t?^,^, 
by u^, and (I2.38P by X^, adding the resulting equations and integrating in time. Developing 
the very same calculations as for Lemma |3.8| we conclude (|4.35|) . whence (|4.28p . as well, for a 
constant independent of ^ > 0. Likewise, arguing by comparison in the equations satisfied by 
1?^ and i}s,fi one obtains that there exists C > such that for all > 

\\9tw^\\L2(^f)^T.v,) + ||t7^||l2(o,t;W') + \\dtZfi\\L'2{o,T:m{rcy) ^ . (4.36) 
Further, we are in the position of proving the following crucial estimate 

\\w^i\\l^(0,T;H) + \\z^i\\L^{0,T;L^{rc)) ^ (4.37) 

for a constant independent of /i > 0. Indeed, let us test (I3.13|) . with e = 0, by w;^ = ^^{i)^), 
(|3.14p . with e = 0, by = ^^(i?^^^), add the resulting relations and integrate on some time 
interval (0,t), t € (0, T] (note that these estimates may be performed rigorously since € 
L^{0,T;V) and G L'^{0,T; H^iVc)) for all ^ > 0). Easy calculations lead to 



(4.38) 



l\\Mt)\\H+ [ / VT?^V^«^ + ^||z^(t)||i2(r,)+ / / Vi?,,^Vz^ 
where 

^14 = ^ X ' ^^■^('^*"'^)ll^^'l - '^II"a'II/^i(o,T;W) + ^ ^ II'^/^IIh ' (4-39) 

h5= f I Mw^,\< I \\h\\H\\w^,\\H, (4.40) 
JO in JO 

-^16 = / / \dt\{X^,)\\z^,\ <C ||9fX^||i2(r^) (llX^llioo(r^) + l) ||2/,||L2(r^) 

JO JTc JO 

< C'll'^/ill/fi(o,T;L2(rc)) + C" / (ll^Mlli°°(rc) + ^ 





(4.41) 



I^MllL2(r,) 



the first inequality in (j4.4ip ensuing from (|2.12p . Now, we remark that the second and the fourth 
term on the left-hand side of (14.380 are nonnegative thanks to ()4.24p and the monotonicity of 
Combining the fact that k takes positive values (cf. (|2.H6p ) with the latter monotonicity 
argument (indeed, it can be easily checked that the trace Wfj^lr^ of on fulfils w^lr^ = 



Bonetti, Bonfanti, Rossi / Thermal effects in contact with adhesion 



33 



-^(^mItc))! conclude that the fifth term (in the l.h.s. of (j4.38p ) is nonnegative as well. Thus, 
we collect (jT^ - ([OT]) : recahing (|2.H8p . estimates (HH]) and (fO|) on the initial data w^^ and 

as well as estimates (I05]l and (lOHl) (the latter yields a bound for in L^{0,T; L'^iTc))), 
and applying the Gronwall lemma, we end up with (|4.37|) as desired. 

All of the above estimates, the Ascoli-Arzela theorem, |28| Them. 4, Cor. 5], and standard 
weak compactness results yield that there exist a subsequence of {{w^,'&^, Zf^, '&s,fi, Ufi,^fi,(,^i, 
(which we do not relabel) and functions {w,'d, z,'ds,u,X, ^,r]) for which the first of (I4.13p . the 
first of (fiTi|) . convergences (|il9]) - (fi:22]) and 

w^^*w mL°^{0,T;H)nH\0,T-V'), 
w^^w in C\[^,T]-V') , 

z^^*z mL^{0,T;L\T,))nH\0,T;H\T,y), 

n 1 , (4-43) 

z^^z mC\[0,T];H\r,y) 

hold as /i \ 0. Arguing in the very same way as in the proof of Proposition 14.41 we find that the 
eight-uple {w,i}, z,i}s,u,X,C,r]) satisfies equations (IT^ . (ITMl) . (IZHBIl - (1X101) . Furthermore, 
combining (jT3l) . ([i^ . (I4.25|) . (|4.42p and (I4.43|) . we conclude that the quadruple {w,z,u,X) 
complies with the initial conditions (|2.28p - (|2.3ip . 

Moreover, as already pointed out in [131 Sec. 4] (again on the basis of [21 Lemma 1.3, p. 42]), 
to conclude (j2.33p it is sufficient to prove 

limsup / / w^i?^ < / -wd (4.44) 
fi\o Jo Jn Jo Jn 

(analogously for proving (|2.35p ). This can be shown by arguing in the very same way as for (j4.34p . 
so we refer the reader to the calculations developed in the proof of Proposition 14.41 

Finally, we recall that estimate (j4.3ip holds for a constant independent of /.i > 0. Further, 
we note that, by definition of .if^ and convergences (I4.13P and (I4.42p . there holds 

l^J.{w^,) = i}^ - fj.Wf_, ^ i} in L^(0, r;i?) , 

so that, by the definition (13. 2p of 7^, — Pfj,{w^) — > in L^(0, T; H). Hence, in view of (I4.42p 

p^{w^,)^w mL\0,T;H). (4.45) 

Using that, by (|3.5p . 3*^{wp) > j* {pfj_{wi^i)) a.e. in 17, we conclude that for all to S (0,T) and 
r > such that (to — to + r) C (0, T) there holds 



/ j*{w) < liminf / j*{p^{w^)) < liminf / fAw^j) < 2rC , 



where again the first inequality follows from ()4.45p and weak lower semicontinuity of the integral 
functional induced by j* . With the same Lebesgue point argument used for proving (j4.3ip . we 
infer (|2.2Up (a completely analogous argument yields (|2.23p ). In the end, we point out that (|2.2Up 
f (12:23]) . respectively), (12^33]) resp.), and (l2Tl2]l yield an estimate for 'd in L°°(0, T; L'^{n)) 

(for in L~(0,r;Li(rc)), resp.). n 
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A Appendix 

The following result shows how the coercivity property ()2.H2p translates in terms of the Yosida 
approximation of the functional j* . 

Lemma A.l. Assume (l2JlT]) - (f2li2]l . Then, 

1. there exists a constant C2 > (Ci being the same constant as in (i2Jl2]l ) such that 

yfi>0,ueV: fi\\^^{u)\\jj+j;{^^{u))>Ci\\u\\mn)-C2; (A.l) 

2. there exists a constant C2 > (Ci being the same constant as in (f2Jl2]) ; such that 

Va.>0, vgH\T,) : ^||if^(^;)||i2(r,)+j;(i?'MW)>^^ill^llLi(r.)-^;. (A.2) 

Proof. We shall just prove (|A.ip . the proof of ()A.2p being completely analogous. For a given 
u ^ V, let us put w := ^^{u): it follows from the definition of and from ()3.2p that 
u — jiw € 7(p^(w)), whence 

p^{w) e(.{u- iiw) . (A.3) 

Therefore, one has 

jliw) > j* {p,,{w)) > Ci||n - pw\\Li(n) - C2 

> C'i||^*IIli{!^) - I^M^^^WMIh - C2 (A.4) 

> Ci\\u\\li^q) - nWwWjj - ^C^ - C2 , 

where the first inequality follows from (13. 5p . the second one from (l2.H2p . the third one from the 
Holder inequality, and the last one from trivial computations. Hence, we conclude (jA.ip . n 

We may now give the 

Proof of Lemma 14.11 We shall just develop the construction of the sequences {w^} and {"i?^}, 
the proof of the second part of the statement being completely analogous. For every > 0, we 
define € ^ as the solution of the variational equation 

/ w^v + p \7w^Vv = / wqv for all v G V. 
Jn Jn Jn 

Arguing in the very same way as in the proof of [3l Lemma 2.4], we find that (|4.iP " (|4.3p hold. 
By (jni) . we have 

which ensures that i!}^ is in V as well, since 7^, is Lipschitz continuous. Then, (j4.5p immediately 
ensues. □ 
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